Measurements of Cosmic-ray Low-energy Antiproton and Proton Spectra in a 
Transient Period of the Solar Field Reversal 
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The energy spectra of cosmic-ray low-energy antiprotons (p's) and protons (p's) have been measured 
by BESS in 1999 and 2000, during a period covering the solar magnetic field reversal. Based on these 
measurements, a sudden increase of the p/p flux ratio following the solar magnetic field reversal was 
observed, and it generally agrees with a drift model of the solar modulation. 

PACS numbers: 98.70.Sa, 96.40.Kk, 95.85.Ry 



Much of the real underlying physics of the Sun is re- 
lated to the 22 yr solar magnetic cycle with recurrent 
positive and negative phases. The magnetic field polar- 
ity reverses when the solar activity is maximum, and the 
global magnetic field profile reverses in the heliosphere. 
The most recent field reversal happened in the beginning 
of 2000. The solar modulation of cosmic rays is caused 
by expanding solar wind, which spreads out locally irreg- 
ular magnetic field and therefore modifies energy spectra 
of the cosmic rays entering the heliosphere. The posi- 
tive and negative particles drift in opposite directions, 
taking different routes to arrive at the Earth in the he- 
liospheric magnetic field. The charge-sign dependence is, 
therefore, a natural consequence jjj], on top of the com- 
mon time dependent change in the overall modulation. 
It explains alternate appearances of "flat" and "peaked" 
periods in neutron monitor data around solar minima. 
Recent works ||-||] a l so indicated that the drift produces 
non-negligible effects between the positive and negative 
particles even during the high solar activity. This view 
is supported by measurements of temporal variation in 
cosmic-ray ratios ||, such as electrons to helium nuclei 
(He) and electrons to protons (p's), where the largest 
variation is associated with the solar magnetic field re- 
versal. Among various cosmic-ray pairs, antiprotons (p's) 
and p's are ideal || for understanding drift effects under 
the change in overall modulation level, because they are 
different only in charge sign. 

In the last solar minimum period, the BESS experi- 
ment revealed that the cosmic-ray p spectrum has a dis- 
tinct peak around 2 GeV ||], which is a characteristic 
feature of secondary p's produced by cosmic-ray interac- 
tions with interstellar (IS) gas. It has become evident 



that p's are predominantly secondary in origin, because 
several recent calculations of the secondary spectrum ba- 
sically agree with observations in their absolute values 
and spectral shapes j|PM^|. 

We report here new measurements of cosmic-ray p and 
p fluxes and their ratios in the energy range from 0.18 
to 4.2 GeV collected in two BESS balloon flights carried 
out in 1999 and 2000, when the solar activity was max- 
imum. Based on the solar magnetic field data ||, the 
Sun's polarity reversed between these two flights ]To| |. 
With our full set of data PJll"|-|l3|], we discuss the tem- 
poral variation of the p flux and p/p ratio covering the 
solar minimum, the maximum, and the field reversal. 

The BESS spectrometer was designed fl4| , [[5[ and de- 
veloped p^l9|] as a high-resolution spectrometer. A uni- 
form field of 1 Tesla is produced by a thin superconduct- 
ing solenoid |2(| , and the field region is filled with track- 
ing detectors. This geometry results in an acceptance of 
0.3 m 2 sr. Tracking is performed by fitting up to 28 hit- 
points in drift chambers, resulting in a magnetic-rigidity 
(R = Pc/Ze) resolution of 0.5% at 1 GV. The upper 
and lower scintillator-hodoscopes jl9| provide time-of- 
flight and two dE/dx measurements. Time resolution 
of each counter is 55 ps, resulting in a 1//3 resolution 
of 0.014, where (3 is defined as particle velocity normal- 
ized by the speed of light. The instrument also incor- 
porates a threshold-type Cherenkov counter jll| with a 
silica-aerogel radiator (n=1.02) that can reject e~ and 
/i - backgrounds by a factor of 4000 and identify p's from 
such backgrounds up to 4.2 GeV. In addition to biased 
trigger modes |R]Jlj]] selecting negatively-charged parti- 
cles preferentially, one of every 60 (30) first-level trig- 
gered events were recorded to provide unbiased samples 
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FIG. 1. The identification plots of p events for (a) 1999 
and (b) 2000 flights. The dotted curves define the p mass 
bands. 

in 1999 (2000). Note that the first-level trigger was pro- 
vided by a coincidence between the top and bottom scin- 
tillators, with the threshold set at 1/3 of the pulse height 
from minimum ionizing particles. 

The experiments were carried out in northern Canada, 
Lynn Lake to Peace River, where the geomagnetic cutoff 
rigidity ranges from 0.3 to 0.5 GV. Data for flux mea- 
surements were taken for live-time periods of 100,471 and 
90,765 sec at altitudes above 34 km (residual air of 4.3 
and 5.0 g/cm 2 on average) in 1999 and 2000, respec- 
tively. Data during ascent were also collected at various 
altitudes, in order to estimate a background from atmo- 
spheric secondary p's. 

Analysis was performed in the same way as described 
in Ref. We applied the same selection criteria for 

p's and p's because non-interacting p's behave like p's 
except for deflection in the symmetrical configuration of 
BESS. Figure [l] shows versus R plots for surviving 
events (for 1999, only the negative rigidity side is shown). 
We see clean narrow bands of 668 and 558 p candidates 
at the exact mirror position of the p's for the 1999 and 



2000 data, respectively. The p bands are slightly con- 
taminated by e~ and backgrounds due to inefficiency 
of the aerogel Cherenkov counter, while contamination 
on the p band was negligible. We estimated the fraction 
of e" and fj,- backgrounds to be %(0 %), 0.6 % (0.5 
%), and 2.3 % (1.5 %), respectively, at 0.3, 2, and 4 GeV 
for 1999 (2000) data. Other backgrounds such as albedo, 
mis-measured positive-rigidity particles, and "re-entrant 
albedo" were found to be negligible. 

The survival probabilities for p's and p's to traverse 
residual air and the instrument without interactions and 
to pass through all the selection criteria described above 
were estimated by a Monte Carlo (MC) simulation. The 
MC code was tuned and verified by comparing the simu- 
lation with an accelerator beam test of the BESS detec- 
tor |^lf . The systematic uncertainty of the p interaction 
losses were reduced to 5 %. Estimates of the background 
from atmospheric secondary p's were calculated |22]-|2^] 
with the p and He fluxes measured by BESS as input. 
The subtraction amounts to 15 ± 5 % (17±6 %), 23 ± 
2 % (30 ± 2 %), and 26 ± 8 % (28 ± 9 %) for 0.3, 2, 
and 4 GeV in 1999 (2000), respectively, where the errors 
correspond to the maximum difference among the cal- 
culations. The background from atmospheric secondary 
p's was estimated according to the calculation p5j . The 
parameters used in the calculation were tuned and ver- 
ified |2(| by comparing the calculated fluxes at various 
air depths with the atmospheric growth curve of p's mea- 
sured by BESS. The subtraction amounts to 27 ± 4 % (55 
± 8 %), 10 ± 2 % (22 ± 3 %), and 4 ± 1 % (8 ± 1 %) for 
0.2, 0.5, and 1 GeV in 1999 (2000), respectively, where 
the errors correspond to the maximum difference between 
observed and calculated atmospheric growth curve of the 
p's. Table | gives the resultant fluxes of p's and p's and 
their flux ratios, p/p, at the top of the atmosphere (TOA) 
with the statistical (first) and systematic (second) errors. 

Figure § shows the p and p fluxes in 1999 and 2000, 
together with the data || taken in 1997 at the solar min- 
imum period. The error bars represent quadratic sums of 



TABLE I. Antiproton fluxes (in 10 2 m 2 s x sr 1 GeV 1 ), proton fluxes (in 10 2 m 2 s 1 sr x GeV x ), and p/p ratio (in unit 
of 10" 5 ) at TOA. T (in GeV) defines the kinetic energy bins. Np is the number of observed antiprotons. 
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FIG. 2. The (a) p and (b) p fluxes at TOA measured by BESS in 1999 and 2000 together with the previous data in 1997. The 
curves represent the calculations by the drift model Q at a period of 10°(+), 70°(+), and 70° (— ), which respectively represent 
solar minimum at positive phase (dashed line), solar maximum at positive phase (dash-dotted line), and solar maximum at 
negative phase (solid line). The dotted curves represent the calculations by the spherically symmetric model Effl. 



the statistical and systematic errors, as well as in the 
following figures. The dotted curves represent the modu- 
lated spectra according to the standard spherically sym- 
metric approach 27 1, in which the modulation is char- 
acterized by a single parameter (cf>) irrespective of the 
Sun's polarity. In each year, <fi was determined to fit the 
p spectrum measured by BESS assuming the IS spectrum 
in Ref. (3). The same (f> was applied to modulate the p 
flux. In both plots, dashed, dash-dotted, and solid curves 
represent modulated spectra according to a steady-state 
drift model || in which the modulation is characterized 
by a tilt angle (a) of the heliospheric current sheet and 
the Sun's magnetic polarity (denoted as +/—). Note that 
periods at 10°(+), 70°(+), and 70°(-) roughly corre- 
spond to the years of 1997, 1999, and 2000, respectively. 
However, it is shown in Fig. ^ that the actual level of 
overall modulation at the date of observation was lower 
(corresponding to a smaller value of a) in 1999 and higher 
(larger a) in 2000. Thus, the predicted effect of the drifts 
in Fig. ||(a) might be hidden under the change in overall 
modulation level. The drift model needs to be fine-tuned 
to reproduce the p and p spectra. The retardation of the 
modulation effect due to the gradually expanding helio- 
spheric current sheet p8[ might also be concerned. 

Figure^ shows the p/p ratios in 1999 and 2000 in which 
the above mentioned difficulty is partly canceled. In gen- 
eral, p's suffer less modulation than p's because of their 
relatively hard IS spectrum. Thus, the spherically sym- 
metric model predicts an increase in p/p towards the so- 
lar maximum period p9fl , mainly due to a decrease of 
the p flux. The drift model predicts a much larger in- 
crease in p/p reflecting the charge-sign dependence. It 
should be noted that a recent independent work ||] has 
shown qualitatively the same feature in p/p while these 
two drift model calculations were different in IS spectra 



and some modulation parameters. As shown in Fig. |^, 
the measured p/p ratios are in better agreement with the 
drift model than the spherically symmetric model. 

Combining all the previous BESS measurements since 
1993 P,pT|-p^[, annual variations of the p/p ratios (rela- 
tive to the ratio in IS) at energies of 0.2 - 0.4 GeV (closed 
square, denoted as 0.3 GeV), 0.8 - 1.3 GeV (open cir- 
cle, 1.0 GeV), and 1.5 - 2.5 GeV (closed star, 1.9 GeV) 
are shown in Fig. §(a). These data show no distinctive 
year-to-year variation before the solar field reversal indi- 
cated by the shaded area in Fig |(a). The curves repre- 
sent calculations of the p/p time variation according to 
the drift model (dashed) and the spherically symmetric 
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FIG. 3. The p/p flux ratios measured by BESS in 1999 
and 2000 with the previous data in 1997. The curves are 
calculations of the ratio at various solar activity [p|,p7|. 
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FIG. 4. (a) Temporal variation of the p/p ratios (relative to 

the ratio in IS) at 0.3 (closed squares), 1.0 (open circles), and 

1.9 GeV (closed stars) measured by BESS since 1993. Curves 

represent the calculated temporal variation of the ratios (see 

text), (b) Temporal variation of Ncl and a. In both plots, 

shaded area indicates the period of the solar field reversal. 

model (dotted) for 0.3, 1.0, and 1.9 GeV from top to 
bottom. We obtained those </>'s at the dates of the BESS 
flights (0bess) and then we evaluated tfi's at any mo- 
ments using a linear relation between i^bess and A*cl, 
where Ncl is the monthly averaged count rates of the 
Climax Neutron Monitor |3(j. The historical trend of 
iVcL is indicated by the solid histogram in Fig. ||(b). In 
Ref. H , the p/p ratios are given as a function of the tilt 
angle, a, rather than the actual time. In order to draw 
the dashed curve in Fig. |^(a), we take a as the mean po- 
sition of the maximum latitudinal extent of the current 
sheet |n] taking no account of the retardation of the 
modulation effect. Note that a is allowed to reach 90° 
in the shaded area to avoid discontinuity in p/p at the 
field reversal. The drift model predicts that the p/p ratio 
is stable during the positive phase and shows a sudden 
increase after the positive-to-negative solar field reversal. 
Our measurements generally agree with the drift model. 

As a conclusion, we have measured the temporal vari- 
ation of the p and p fluxes covering the solar minimum, 
the maximum, and the solar magnetic field reversal. We 
observed stable p/p ratios in the positive polarity phase 
through 1999 and a sudden increase in 2000 following the 
solar field reversal. The continuous observation has en- 
abled a crucial test of drift effects in the solar modulation 
including a solar maximum period. It will motivate fur- 
ther development of models with more realistic parame- 
ters jl| and time dependence [jj2|. Furthermore, detailed 
understanding of solar modulation is inevitably impor- 
tant for investigating the origin of p's |33] , especially in 



the very low energy region below 0.5 GeV, where slightly 
excessive p fluxes relative to the theoretical calculations 
were observed during the solar minimum 
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